Astronomy & Astrophysics manuscript no. 16671 


©ESO 2013 


January 12, 2013 





WASP-39b: a highly inflated Saturn-mass planet orbiting a late 

G-type star 

F. Faedi^ S. C. C. Barros', D. R. Anderson^, D. J. A. Brown^, A. Collier Cameron^, D. PolIacco\ I. Boisse"^-^, G. 

Hebrard^ ^, M. Lendl^ T. A. Lister^ B. Smalley^, R. A. Street^ A. H. M. J. Triaud^ J. Bento^ O. W. Butters'", B. 
Enoch^ F Bouchy4'6, C. A. Haswell'\ C. Hellier^, F R Keenani, G. R. M. Miller^, V. Moulds^ C. Moutoui^, A. J. 
Norton'!, q QueIoz^ A. Santerne^'^, E. K. Simpson^ I. Skillen'^ A. M. S. Smith^, S. Udry^ C. A. Watson^ R. G. 

West'", and R J. Wheatley"^ 

' Astrophysics Research Centre, School of Mathematics and Physics, Queen's University Belfast, University Road, Belfast BT7 
INN. 

e-mail: f . f aedi@qub .ac.uk 
^ Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK 

School of Physics and Astronomy, University of St Andrews, St Andrews, Fife KY16 9SS, UK 

Institut d'Astrophysique de Paris, UMR7095 CNRS, Universite Pierre & Marie Curie, France 
' Centro de Astroffsica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 

Observatoire de Haute-Provence, CNRS/OAMP, 04870 St Michel I'Observatoire, France 
^ Observatoire astronomique de 1' Universite de Geneve, 51 ch. des Maillettes, 1290 Sauvemy, Switzerland 

Las Cumbres Observatory Global Telescope Network, 6740 Cortona Drive Suite 102, CA 931 17, USA 
' Department of Physics, University of Warwick, Coventiy CV4 7AL, UK 

Department of Physics and Astronomy, University of Leicester, Leicester, LEI 7RH 
" Department of Physics and Astronomy, The Open University, Milton Keynes, MK7 6AA, UK 

Laboratoire d'Astrophysique de Marseille, 38 me Frederic Joliot-Curie, 13388 Marseille cedex 13, France 

Isaac Newton Group of Telescopes, Apartado de Correos 321, E-38700 Santa Cruz de Palma, Spain 
Received; accepted 

ABSTRACT 

We present the discovery of WASP-39b, a highly inflated transiting Saturn-mass planet orbiting a late G-type dwarf star with a 
period of 4.055259 ± 0.000008 d. Transit Epoch To = 2455342.9688 ± 0.0002 (HID), of duration 0. 1 168 ± 0.0008 d. A combined 
analysis of the WASP photometry, high-precision follow-up transit photometry, and radial velocities yield a planetary mass of Mpi = 
0.28 ± 0.03 M] and a radius of Rpi = 1.27 ± 0.04 resulting in a mean density of 0.14 ± 0.02 pj. The stellar parameters are mass 
M^, = 0.93 ± 0.03 M0, radius R^, = 0.895 ± 0.23 i?o, and age 9!^Gyr Only WASP-17b and WASP-31b have lower densities than 
WASP-39b, although they are slightly more massive and highly irradiated planets. From our spectral analysis, the metallicity of 
WASP-39 is measured to be [Fe/H] = -0.12 ± 0.1 dex, and we find the planet to have an equilibrium temperature of 1116;^32K. 
Both values strengthen the observed empirical correlation between these parameters and the planetary radius for the known transiting 
Saturn-mass planets. * 

Key words, planetary systems - stars: individual: (WASP-39, GSC 04980-00761) - techniques: radial velocity, photometry 



1. Introduction 

The importance of transiting extrasolar pla nets is related to 
their geometrical configuration (Sackett 1999). Transit geometry 
severely constrains the orbital inclination of the planet, allowing 

* Spectroscopic and photometric data are available in electronic form 
at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via|http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A-l-A/, 



accurate measurements of its mass and radius to be derived. The 
inferred planet's density provides information on the system's 
bulk physical properties, and thus is a fundamental parameter for 
constraining theore tical models of planetary formation , structure 
and evolution (e.g. lGuillotll2005t iFortnev et all 120071; iLiu et all 
2008). 

To date, more than 100 transiting planets have been dis- 
covered, which show a huge range of diversity in their 
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physical and dynamical properti es. For example, th eir mass 
ranges from ~ 5Mffi (Kepler-lOb, 'Batalha et al. 201 1) to about 
12 Mj (XO-3b, Johns-Krull et al. 2008; and Hebrard et al. 2008\. 
Some planet s have radii that ag r ee with models of irradi- 
ated planets jBurrows et al.l l2007t iFortnev et all l2007l) . while 
others are found to be anom alously large (e.g . WASP- 12b 
Hebb et al.' '2009' and TrES-4b ISouthworthI 1201 oi iTorres et al.l 
2008, Mandushev et al. 2003). The diversity in exoplanet den- 
sities, and hence in their internal compositions, is particularly 
noticeable at sub-Jupiter masses. For example, some exoplanets 
have very high densities and ar e thought to have a rocky/ice core 
(e.g. HD 149026b, ppi ^ Ipj 
such as TrES-4b (ppi - 0.17 pj. 



Sato et al."2005), while systems 
Mandushev et al, 2007). WASP- 



ITbJppi = 0.06pj, lAnderson et al l feoiObt [Anderson et al.l 



l20TT bl). WASP-31b (o„i = 0.132 ot. lAnderson et al.ll2010al) . and 
Kepler-7b (ppi - 0.1 3 pi. [Latham et alJ 2O10l) are examples of 
planets with puzzlingly low densities that challenge s tandard 
evolu t ionary theories in re producing their radii ( Fortnev et alj 
120071 iBurrows et al.ll2007h . To assess the inflation status of a 
system, generally planetary radii are compared to tabulated val- 
ues from models (e .g. Fortnev et al. 2007, Burrows et al. 2007, 
iBaraife et al.ll2008h . However, the radius depends on multiple 
physical properties, such as the stellar age, the irradiation flux, 
the planet's mass and age, the atmospheric composition, the 
presence of heavy elements in the envelope or in the core, the 
atmospheric circulation, and also on any source generating ex- 
tra heating in the planetary interior Different mechanisms have 
been proposed to explain the anomalously large radii, such as 
tidal h eating due to unseen companion s pumping up the eccen- 
tricity (iBodenheimeret al.ll200l[ '2003'). kinetic heating due to 
the breaking of atmospheric waves (Guillot & Showman 20021), 
enhanced atmospheric opacity ( Burrows et al. 2007) and semi- 
convection (Chabrier & Baraffd 12007.) . While each individual 
mechanism would presumably affect all hot Jupiters to some 
degree, they cannot explain th e entirety of the o bserved radii 
tFortnev & Nettelmann 2010, Baraffe et all 1201 Ol) . More com- 
plex thermal evolution models are necessary to fully understand 
their cooling history. For these planets, the dominant source of 
energy is a function of the orbital separation and the spectral 
type of the host star, while its dependency on the planetary mass 
and age is negligi ble (with the excepti on of very young plan- 
ets). More recentlv. batvgin et al.l(l201 ih performed calculations 
of the thermal evolution of gas giants and suggest that the extra 
energy needed to ex plain the radius infl ation comes from stellar 
irradiation (see also lLauehlin et all201 ih . The proposed mecha- 
nism (Ohmic heating), based on the interaction of ionised alkali 
metals in the planetary atmosphere with the planet's magnetic 
field, along with the atmospheric heavy element content, could 
provide a universal explanation of the currently measured radius 
anomalies. 

Although the majority of the known exoplanets are short- 
period, Jupiter-mass planets, more recently an incr easing num- 
ber of Saturn-like planets have been discovered (e.g. lEnoch et alj 
l20Tlh . and have encouraged studies of planetary properties 
and their statistical analysis, searching fo r possible correla- 
tions be tween planetary parameters (e.g. lEnoch et"aLl 1201 It 
lAnderson et al. 2010a; Hartman et al. 20111). To date, 27 tran- 
siting planets have been discovered with masses in the range 
0.15 Mj < Mpi < 0.6 MjB similar to Saturn (Msatum = 0.229 M,, 
IStandish|[T995l) . The detection and characterisation of signifi- 
cantly more bright short-period transiting systems is one of the 
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keys to improving our understanding of planetary structure and 
evolution. 

Here we describe the properties of WASP-39b, a new tran- 
siting Saturn-mass planet discovered by the SuperWASP survey. 
The planet host star WASP-39 belongs to the constellation of 
Virgo and thus resides in an equatorial region of sky, which 
is monitored by the SuperWASP-North and WASP-South tele- 
scopes simultaneously. WASP-39b is the third least dense planet 
(Ppi = 0.14 ± 0.02 pj) discovered from a ground-based transit 
survey, and belongs to the sample of highly inflated gas giant 
planets. It provides observational evidence for the mass-radius 
relation of planetary systems in a poorly sampled region of the 
parameter space. 

We present follow-up observations of the new system which 
establish the planetary nature of the transiting object detected by 
SuperWASP. High precision, high signal-to-noise light-curves 
have been obtained using both the Faulkes Telescope North 
(FTN), and the Euler telescope, and radial velocity measure- 
ments using the SOPHIE (1.93-m OHP) and CORALIE (Swiss 
1.2-m) spectrographs. 

The paper is structured as follows: in § 2 we describe the 
observations, including the WASP discovery data and the photo- 
metric and spectroscopic follow-up. In § 3 we present our results 
for the derived systems parameters, and the stellar and planetary 
properties. Finally, we discuss the implication of the discovery 
ofWASP-39b in §4. 



2. Observations 

ISWASP J142918.42-032640.1 (2MASS 14291840-0326403), 
hereafter WASP-39, has been identified in seve ral northern sky 
catalogues which provide broad-band o ptical (jZacharias et al.l 
120051) and infra-red 2MASS magnitudes (ISkrutskie et al.ll2006l) 
as well as proper motion information. Coordinates, broad-band 
magnitudes, and proper motion of the star are from the NOMAD 
catalogue and are given in Table [1] 



2.1. SuperWASP observations 

The WASP North and South telescopes are located in La Palma 
(ING - Canaries Islands) and Sutherland (SAAO - South Africa), 
respectively. Each telescope consists of 8 Canon 200mm f/1.8 
focal lenses coupled to e2v 2048 x 2048 pixel CCDs, yielding 
a field -of-view of 7.8 x 7.8 square degrees with a pixel scale of 
13.7" (iPollacco et alJl2006l) . 

WASP-39 is a V = 12.11 star located in an equatorial re- 
gion of sky monitored by both WASP instruments, significantly 
increasing the observing coverage on the target. In January 
2009, the SuperWASP-N telescope underwent a system upgrade 
that improved our control over the main sources of red noise, 
such as temperature-dependent focus changes (Barr os et al.l 
This upgrade yielded data of unprecedented high quality, 
and increased the number of planet candidates, flagged in the 
archive, in particular those with longer period and lower mass 
(e.g. WASP-38b, BaiTos et al. 201 1 ; and WASP-39b). 

WASP-39 was routinely observed between 2006 July 1 to 
2010 July 26, with a total of 11 WASP fields and 40531 pho- 
tometric points. Over the five WASP seasons only three points 
were taken in 2006 and none during 2007. However, the same 
field was observed again in 2008, 2009 and 2010 after both 
WASP telescopes began observing an overlapping equatorial re- 
gion. 
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Table 1. Photometric properties of WASP-39. The broad 
band magnitudes and proper motion are obtained from the 
NOMAD 1.0 catalogue. 



1 00 



Parameter 


WASP-39 


RA(J2000) 


14:29:18.42 


Dec(J2000) 


-03:26:40.1 


B 


12.93+0.25 


V 


12.11±0.13 


I 


11.34±0.08 


J 


10.663+0.024 


H 


10.307+0.023 


K 


10.202+0.023 


yUff(mas/yr) 


-12.2 +3.3 


y"i(mas/yr) 


2.8+3.5 
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Fig. 1. Upper panel: Discovery light curve of WASP-39b phase 
folded on the ephemeris given in Table |4] Lower panel: binned 
WASP-39b light curve. Black-solid line, is t he best-fit tran- 
sit model estimated using the formalism from iMandel & Agoil 
(l2002h . 
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Fig. 2. Upper panel: Phase folded radial velocity measure- 
ments of WASP-39 obtained combining data from SOPHIE 
(filled-circles) and CORALIE (open-squares) spectrographs. 
Superimposed is the best-fit model RV curve with parameters 
from Table |4] The centre-of-mass velocity for each data set 
was subtracted from the RVs (ysoPHiE - -58.4826kms"'and 
TcoRALiE = -58.4708 km s"'). Lower panel: we show the bisec- 
tor span measurements as a function of radial velocity, values 
are shifted to a zero-mean (< Yspan ^sophie— 

-0.032kms-i, 

< Vspan >coRALiE- -0.05 1 km s"' ). The bisector span shows 
no significant variation nor correlation with the RVs, suggesting 
that the signal is mainly due to Doppler shifts of the stellar lines 
rather than stellar profile variations due to stellar activity or a 
blended eclipsing binary. 



All data were pro cessed with a custom -built reduction 
pipeline described in iPoUacco et al.l (l2006h . The resulting 
light curves were analysed using our implementation of the 
Box Least-Squares and Sy s Rem detrending alg orithms (see 
Collier Cameron et al.l 120061: iKovacs et al.l 120021: [Tamuz et al., 



2005h . to search for signatures of planetary transits. When com- 



bined, the SuperWASP light curves showed a characteristic peri- 
odic dip in brightness with a period of P = 4.055 days, duration 
Tin ~ 168 minutes, and a depth ~ 25mmag. Figure [1] shows the 
discovery photometry of WASP-39b phase folded on the above 
period, along with the binned phased light curve. A total of 10 
partial or full transits were observed, with an improvement inx^ 
of the box-shaped model over the flat light curve of Ax^ - 825. 
We evalua ted a signal-to-(red) noise value of the data following 
iPont et al.l (I2006X and determined SNred = 1 1.65. 



2.2. Spectroscopic follow-up 

WASP-39 was first observed during our follow-up campaign in 
April 2010 at Observatoire de Haute-Provence (OHP). During 
our program we have obtained follow-up spectroscopy and es- 
tablished the planetary nature of WA SP-39b together w ith three 
addi tional systerns: WA SP-37b (Si mpson et alJl201 ih. WASP - 
38b (iBarros et al.l201lh and WASP-40b (lAnderson et al.l201 lal) . 
Between 2010 April 8 and June 11, we obtained eight radial 
velocity measurements for WASP-39 using SOPHIE, the fiber- 
fed echell e spectrograph mounted on the 1 ■93-m te lescope at 
the OHP dPerruchotet all 120081; iBouchv et al.ll2009h . We used 
SOPHIE in high efficiency mode (R = 40,000) and obtained 
observations with very similar signal-to-noise ratio (~30), to 
minimise systematic errors arising from the known Charg e 
Transfer Inefficiency effect of the CCD dBouchv et all l2009l) . 
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Table 2. Radial velocity (RV) and line bisector span (Vspan ) mea- 
surements of WASP-39. 



BJD 


RV 


o-Rv 


V 

' span 


Instrument 


-2,400,000 


(km s-i) 


(kms-') 


(kms-') 




55321.7484 


-58.455 


0.011 


-0.017 


CORALIE 


55327.7544 


-58.523 


0.022 


-0.047 


CORALIE 


55359.5932 


-58.511 


0.011 


-0.067 


CORALIE 


55361.6552 


-58.502 


0.035 


-0.128 


CORALIE 


55362.5676 


-58.471 


0.013 


-0.050 


CORALIE 


55365.5935 


-58.455 


0.012 


-0.048 


CORALIE 


55376.6394 


-58.534 


0.012 


-0.071 


CORALIE 


55378.6265 


-58.473 


0.009 


-0.036 


CORALIE 


55380.6376 


-58.522 


0.012 


0.006 


CORALIE 



1 00 



50 - 



-50 - 



-1 00 - 



5320 5340 5360 5380 

Time (BJD - 2450000.0) 



55304.4704 


-58.475 


0.013 


0.000 


SOPHIE 


55305.4374 


-58.429 


0.015 


-0.034 


SOPHIE 


55323.4389 


-58.565 


0.017 


-0.031 


SOPHIE 


55329.4184 


-58.472 


0.016 


-0.019 


SOPHIE 


55331.4301 


-58.522 


0.014 


-0.040 


SOPHIE 


55334.4624 


-58.418 


0.017 


-0.036 


SOPHIE 


55336.4477 


-58.503 


0.016 


-0.035 


SOPHIE 


55368.3955 


-58.514 


0.017 


-0.064 


SOPHIE 



although this is now corrected by the data reduction software. 
Wavelength calibration with a Thorium- Argon lamp were per- 
formed every ~ 2 hours, allowing the interpolation of the spec- 
tral d rift of SOPHIE (< 3 m s ' per hour; see iBoisse et al.l 
I2OIOI) . Two 3" diameter optical fibres were used, the first cen- 
tred on the target and the second on the sky to simultane- 
ously measure the background to remove contamination from 
scattered moonlight. During our observations the contribution 
from scattered moonlight was negligible as it was well shifted 
from the target's radial velocity. Nine additional radial veloc- 
ity measurements were obtained using the CORALIE spectro- 
graph mounted on the 1 .2-m Euler S wiss telescope at La Silla, 
Chile dBaranne et al. 1996; Oueloz et al. 2000; Pepe et al.l2002h . 
Observations were obtained with a signal-to-noise of ~30, dur- 
ing grey/dark time to minimise moonlight contamination. The 
data were processed with the SOPHIE and CORALIE standard 
data reduction pipelines, respectively. Radial velocity uncertain- 
ties were evaluated, inc luding known syst ematics such as guid- 
ing and centring errors (iBoisse et al.ll2010]) . and wavelength cal- 
ibration uncertainties. All spectra were single-lined. 

We computed the radial velocities from a weighted cross- 
correlation of each spectr um with a numerical m ask of spectral 
type G 2, as described in iBaranne et al.l (Il996h and iPepe et aP 
(|2002|) . The cross-correlation with masks of different spectral 
types (FO, K5 and M5) produced similar radial velocity varia- 
tion, rejecting a blended eclipsing system of stars with unequal 
masses as a possible cause of the variation. Radial velocity mea- 
surements and line bisector (Vspan) are given in Table|2] and plot- 
ted with the best-fit Keplerian model in Figure |2] SOPHIE data 
are plotted as filled circles and CORALIE data as open squares, 
and both data sets are offset with respect to the radial veloc- 
ity zero point, ysoPHiE and ycoRALiE, respectively (see Table |4|i. 
No significant correlation is observed between the radial veloc- 
ity and the line bisector, suggesting the signal's origin as plan- 
etary rather than due to a blende d eclip sing binary system or to 
stellar activity (see lOueloz et al.ll2001l) . The RMS for SOPHIE 
and CORALIE radial velocity residuals to the best-fit model are 




5320 5340 5360 5380 

Time (BJD-2450000.0) 



Fig. 3. Upper panel: The bisector span measurements as a func- 
tion of time (BJD-2450000.0), Vspan values are shifted to a zero- 
mean as in Figure|2] Lower panel: Residuals from the RV orbital 
fit plotted against time. 



RMS SOPHIE = 26 m s"' and RMS coralie = 16 m s"\ with the 
higher RMS value of the SOPHIE measurements dominated by 
two discrepant measurements, RVi = 48 m s"' and RV2 = 43 m 
s"'. Eliminating the most discrepant of these, the RMS sophie 
becomes 19m s comparable to that of CORALIE. We also 
plot in Figure |3] the bisector span measurements as a function 
of time; and lower panel, the residuals from the RV orbital 
fit against time. The bisector span values are shifted to a zero 

mean (< Vspan >S0PHIE= -0.032kmS"', < Vspan >coralie = 

-0.05 1 km s"'), to better compare the two data sets. 
2.3. Photometric follow-up 

To allow more accurate light curve modelling and thus refine 
the photometric parameters, we obtained two high signal-to- 
noise transit light curves of WASP-39b. All photometric data 
presented here are available from the NStED database 0. The 
first full transit was observed with the LCOG10 2m Faulkes 
Telescope North (FTN) on Haleakala, Maui Hawai'i, on the 
night of 2010 May 18. The fs03 Spectral Instruments camera 
was used with a 2 x2 binning mode giving a field-of-view of 1 0' x 
10' and a pixel scale of 0.303"/pixel. Data were taken through 



- http://nsted.ipac.caltech.edu 
' http://lcogt.net 
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Orbital phase 

Fig. 4. FTN z-band and Euler r-band follow-up high signal- 
to-noise photometry of WASP-39b during the transit. The Euler 
light curve has been ofsetted from zero by an arbitrary amount 
for clarity. The data are phase-folded on the ephemeris from 
Table |4] Superposed (black-solid line) is the be st-fit transit 
model estimated using the formalism from Ma ndel & Agoll 
(l200l . Residuals from the fit are displayed underneath. 



a Pan-STARRS-z filter and the telescope was defocussed dur- 
ing observations to prevent saturation and allow longer expo- 
sure times to be used. T he observations were pre-processed us- 
ing the WASP Pipeline dPollacco et alj|2006[) to perform master 
bias and flat construction, debiassing and flatfielding. Aperture 
photome try was performed with the the DAOPHOT package 
(IStetsonlll987l) within the IRAFQ environment using a 1 3 pixel 
radius aperture. The differential photometry was performed rel- 
ative to 7 comparison stars that were within the FTN field of 
view. Additional high signal-to-noise photometry was obtained 
in the Gunnr filter with the Euler-Swiss telescope on 2010 July 
9, only covering a partial transit. Conditions were variable, with 
seeing ranging from 0.6 " to 1.7 ". The Euler telescope employs 
an absolute tracking system, which keeps the star on the same 
pixel during the observation by matching the point sources in 
each image with a catalogue, and adjusting the telescope point- 
ing between exposures to compensate for drifts. Observations 
were obtained with a slightly defocused (0. 1 mm) telescope. All 
images were corrected for bias and flat field eff'ects and the light 
curve was obtained by performing relative aperture photometry 
of the target and one bright reference star. 

Both the FTN and Euler light curves are shown in Figure 2] 
The photometry confirms the transit which phases with the 
ephemeris derived from the WASP discovery photometry. 



* IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Table 3. Stellar parameters of WASP-39 from Spectroscopic 
Analysis. 



Parameter 


Value 




5400 ± 150 K 


log,? 


4.4 ± 0.2 




0.9 ±0.2 km S-' 


V sin i* 


1.4 ±0.6 km S-' 


[Fe/H] 


-0.12 ±0.10 


[Na/H] 


-0.04 ±0.10 


[Mg/H] 


0.06 ±0.11 


[Al/H] 


0.01 ± 0.08 


[Si/H] 


0.04 ± 0.08 


[Ca/H] 


0.01 ± 0.14 


[Sc/H] 


0.02 ±0.19 


[Ti/H] 


-0.03 ±0.10 


[V/H] 


-0.08 ±0.17 


[Cr/H] 


-0.07 ±0.10 


[Mn/H] 


-0.03 ± 0.20 


[Co/H] 


-0.10 ±0.12 


[Ni/H] 


-0.06 ± 0.09 


log A(Li) 


<0.9 


Mass 


0.93 ± 0.09 Mo 


Radius 


1.00 ± 0.25 Ro 


Sp. Type 


G8 


Distance 


230 ± 80 pc 



Note: Mass and Radius estimate using the lTorres etal.l ( l2010h 
calibration. Spectral Type estimated from T^ff using the table in lGravl 

(20081 . 



3. Results 

3.1. Stellar parameters 

A total of nine individual CORALIE spectra of WASP-39 were 
co-added to produce a single spectrum with a typical S/N ra- 
tio of around 50. The standard pipeline reduction products were 
used in the analysis. To improve the line profile fitting for equiv- 
alent width measurements, the spectrum was smoothed using a 
Gaussian width cr = 0.05 A. For the v sin ;* determination the 
unsmoothed spectrum was used. 

Our ana lysis was performed using the methods given in 
Gill on et alJ m)0%. The H„ line was used to determine the ef- 
fective temperature (Tsff), while the Na i D and Mg i b lines were 
used as surface gravity (logg) diagnostics. The atmospheric pa- 
rameters obtained from the analysis are listed in Table |3] The 
elemental abundances were determined from equivalent width 
measurements of several clean and unblended lines. A value for 
microturbulence (ft) was determined from the Fe i lines using 
the method of Magain (1984). Quoted eiTor estimates include 
those given by the uncertainties in T^s, logg and ft, as well as 
the scatter due to measurement and atomic data uncertainties. 

The projected stellar rotation velocity (vsin/*) was deter- 
mined by fitting the profiles of several unblended Fe i lines. 
A value for macroturbulence (v^ac) of 2.1 ± 0.3 kms~' was 
assumed, based on the tabulation by iGravl (l2008h and an in- 
strumental FWHM of 0.11 + 0.01 A, determined from the tel- 
luric lines around 6300 A. A best-fitting value of vsin/* = 
1.4 ± 0.6 km s"' was obtained. The stellar mass Mj_and_ja- 
dius R j, were estimated using the calibration of iTorres et"all 
(I2010h . 

The non-detection of lithium in the spectrum, the low rota- 
tion rate implied by the v sin /* and lack of stellar activity (shown 
by the absence of Ca ii H and K emission) all indicate that the 
star is relatively old. Unfortunately, the gyrochronological age 
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Table 4. System parameters of WASP-39 
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Fig. 5. Isochrone ti-acks from lDemarque et al"](l2004l) for WASP- 
39 using the best-fit metallicity of [Fe/H] = -0. 12 dex and stellar 
density 1 .297 p©. Solid lines are for isochrones of, left to right: 
1.0,3.0,6.0,9.0, 12.0, 15.0, 18.0, and 20.0 Gyr. Dashed lines ai-e 
for mass ti-acks of, left to right: 1.0, 0.9, 0.8 Mq. 



Parameter (Unit) 
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P{d) 
To (HJD) 
Th (d) " 
Tn = T34 (d) 



AF 
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Ki (m s-') 
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rcoRALiE (km s"') 
e 

(Mo) 
(Rg) 
logg* (cgs) 

P* (Po) 

Mpi (Mj) 
^pi (Rj) 

loggpi (cgs) 

Ppi (Pj) 
a (AU) 



' pl,A=0 



(K) 



4.055259±0.000009 
2455342.9688±0.0002 
0.1168±0.0008 
0.0179+0.0009 

QQ2 1 1+0.0003 
"•"■^^^-0.0004 

441+0.036 

on Qo+0.25 
°' ■0^-0.22 

38±4 
-58.4826±0.0004 
-58.4708+0.0004 
(fixed) 
0.93+0.03 
0.895+0.023 
4.503+0.017 

1 907+0.082 
i-^^ '-0.074 

0.28+0.03 
1.27+0.04 

2 610+" °'" 

0.14+0.02 
0.0486+0.0005 

1116!^ 



" T\a: time between V and 4"' contact 



estimate from the iBame i (l2007l) relation (~ 5!f Gyr) can only 
provide a weak constraint on the age of WASP-39. 

The stellar density of - 1.297^qq^^pj obtained from 
the MCMC analysis was used together with the stellar tem- 
perature and metallicity values, derived from spectroscopy, in 
an interpolation of th e Yonsei-Yale stellar evolution tracks 
(iDemarque et al.ll2004l) . as shown in Figure |5] Using the best-fit 
metallicity of [Fe/H] - -0.12, we obtain a mass for WASP-39 of 
0.86±0.05 Mq and a stellar age of 9^^ Gyr, in agreement with the 
gyrochronological age and a more accurate estimate. The mass 
obtained from the YY-isochrone fit is somewhat less than the 
value from the lTorres et al.l (12010 ^ calibration (as was also found 
in the analysis of WASP-37 stellar parameters; ISimpson et al.l 
1201 Ih . but their 1-cr errors overlap. 

3.2. Planetary parameters 

The planetary properties were determined using a simultane- 
ous Markov-Chain Monte Carlo (MCMC) analysis including 
the WASP photometry, the follow-up FTN and Euler photom- 
etry, together with SOPHIE and CORALIE radial velocity mea- 
surements. A detailed descri pti on of the method is given in 
ICollier Cameron et al.l(l2007l) and lPoUacco et al.l(l2008l) . The pa- 
rameters we used in the fit are: the epoch of mid-transit To, the 
orbital period P, the fractional change of flux proportional to the 
ratio of stellar to planet surface areas AF = R^JR-^, the transit 
duration Tu, the impact parameter b, the radial velocity semi- 
amplitude Ki, the stellar effective temperature T^s, metallicity 
[Fe/H], the Lagrangian elements V^cos u and yfe sin cj (where 
e is the eccentricity and u the longitude of periastron), and the 
systematic off'set velocity y. In this particular case we fitted the 2 



systematic velocities jsophie and jcoralie to allow for instru- 
mental offsets between the two datasets. 

Four diff'erent sets of solutions were considered: with and 
without the main-sequence mass-radius constraint in the case 
of circular orbits and orbits with floating eccentricity. For each 
solution we have included a linear trend in the systemic veloc- 
ity, as a free parameter. Ho wever, we find no s ignificant vari- 
ation. We used the model of lClaretl (l2000l 12004) for the limb- 
darkening in the r-band, for both WASP and Euler photom- 
etry, and in the z-band for FTN photometry. Due to the low 
mass of WASP-39b, the radial velocity data do not off'er con- 
vincing evid ence for an eccentric o rbit. We performed a Lucy 
& Sweeney (iLucv & Sweenevll 19711 Eq. 27) F-test, which indi- 
cates that there is a 54% probabiUty that the improvement in the 
fit produced by the best-fitting eccentricity could have arisen by 
chance if the orbit were truly circular. Moreover, we find that 
imposing the main-sequence constraint has little effect on the 
MCMC global solution. Thus, we decided to adopt no main- 
sequence prior and circular orbit. 

From the above parameters, we calculate the mass M, ra- 
dius R, density p, and surface gravity \ogg of the star (which 
we denote with subscript ★) and the planet (which we denote 
with subscript pi), as well as the equilibrium temperature of the 
planet assuming it to be a black-body 7'pi,A=o and that energy is 
efficiently redistributed from the planet's day-side to its night- 
side. We also calculate the transit ingress/egress times T12/TJ4, 
and the orbital semi-major axis a. These calculated values and 
their 1-cr uncertainties from our MCMC analysis are presented 
in Table |4] The corresponding best-fitting transit light curves are 
shown in Figures [U Figure |4] and the best-fitting RV curve in 
Figure |2] 
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4. Discussion 



We report the discovery of a new transiting extrasolar planet, 
WASP-39b. A simultaneous fit to transit photometry and ra- 
dial velocity measurements gives a planetary mass of 0.28 ± 
0.03 Mj and a radius of 1 .27 + 0.040 Rj which yields a planetary 
density of 0.141 + 0.02 pj. Thus, WASP-39b is the third least 
dense planet identified by a ground-based transit survey. Only 
WASP-17b (jAnderson et al. foiObjp wn ^0 06pi), and WASP- 
31b {pw3i — 0. 132pT. [Anderson etal.l|2010ah have lower densi- 
ties. However, they are sUghtly more massive planets (Mwn - 
0.49 Mj and Mw3i = 0.48 Mj), but highly irradiated, wit h larger 
and ho tter host stars (Anderson et al. 2010b; Anderso n et al.l 
1201 Oah . This implies higher planet equilibrium temperatures for 
WASP-17b and WASP-31b compai-ed to WASP-39b. We find 
WASP-39b to have a highly inflated radius (/?pi = 1 .27 Rj), more 
than 20% larger tha n the j^pi obtained by comparison with the 
iFortnev et M ilQQlh and the iB araffe etal., (2008) models for a 
coreless planet of a similar mass, orbital distanc e and s tellar age. 
For example, tables presented in iFortnev et al.l (12007 ^ predict a 
maximum radius of ~ 1.05 Rj for a 0.24 Mj planet orbiting at 
0.045 AU from a 4.5 Gyr solar-type star. In addition, WASP- 
39 is smaller, cooler and probably older than the Sun. Thus, a 
radius of 1.27 Rj is clearly too large for these models. The fact 
that we do not detect any eccentricity, and that the age of WASP- 
39b's host star is ~ 9 Gyr, suggests that it is unlikely that recent 
tidal circularisation an d dissipation could b e a cause of the large 
radius of WASP-39b (iLeconte et al.1 l20Tol lHansenll20Toh . The 
low metalhcity ([Fe/H] = -0.12 ± 0.1) of the WASP-39b host 
star supports the expected low planetary core-mass. However, 
this will only marginally explain the highly inflated radius of 
WASP-39b. Hence, this l eads to the hypothes is that some addi- 
tional physics is at play (IFortnev et al.ll2009h . An interior heat 
source, replacing the radiated heat from gravitational contrac- 
tion, is needed for the planet to reach thermal equilibrium with 
a larger radi us than theoretically expected. The Ohmic heating 
mechanism ( Batvgin et alj|20lll) . based on the electro-magnetic 
interaction betweent atmospheric winds and the planet's mag- 
netic field, is able to explain the anomalous sizes of close-in 
transiting planets. WASP-39b has a low equilibrium tempera- 
ture (rpi_A=o - 1116 K) and thus ap pears to belong to the 'pL' 
class of planets from IFortnev et al] (12008). We note that, with 
the measured planetary parameters (mass, radius, Teg), WASP- 
39b ag rees with the coreless models in Figure 4 of lBatvgin et al.l 



I 



(1201 

Of the known transiting systems, WASP-39b joins an in- 
creasing number of recently discovered exoplanets with Saturn- 
like masses. With an increasingly large sample of well- 
characterised systems, we can begin to make statistical infer- 
ences as to the ph ysical reasons behind their diverse nature. 
lEnoch et"al] (1201 I h and Anderson et al. (2011 in prep.), show 
that the radii of known low-mass (0. 15-0.6 Mj) giant planets 
strongly correlate with equilibrium temperature and host- star 
metallicity. We have investigated how WASP-39b relates to the 
rest of the sample of Saturn-mass planets. W e calculated the ra- 
dius anomaly "R, as in Laughlin et al. (201 1), for all Saturn-mass 
planets included in the list from Enoch et al. (201 1), and added 
the latest discoveries (for an updated list see the extrasolar planet 
encyclopedia http://exoplanet.eu/). In Figure|6]we plot K against 
the planetary equilibrium temperature (upper panel), and as a 
function of the host star metallicity (lower panel). WASP-39b is 
indicated with a filled circle. Figure |6] shows that K correlates 
wit h the equilibrium tem perature of the planet, as also found 
by iLaughhn et al ] (120111) . For transiting gas giant planets the 
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Fig. 6. Upper panel: The radius anomaly H - Robs - Rpred, ver- 
sus equilibrium temperature for the known Saturn-mass planets. 
Lower panel: % as function of the stellar metallicity [Fe/H] in 
dex. WASP-39b is indicated with a filled circle. 

dominant source of energy is a function of the orbital separa- 
tion and the spectral type of the host star, and depends negligi- 
bly on the planetary mass and age (with the exception of very 
young planets). Thus, the observed "R- Tgq trend supports the re- 
latio n between the inflati on mechanism and the stellar irradiation 
flux (Laughhn et al."2011jD. The masses of planets in the Saturn- 
mass range also appear to correlate with their host star metal- 
hcity dGuillot et al.l 120061; iBurrows et al.H2007ft . such that low- 
density planets are found to orbit su b-solar metallicity st ars (for 
example WASP-21 [Fe/H] = -0.4. iBouchv et al.1 1201 Ol) . while 
higher density planets orbit stars with super-solar metallicit 
(for example WASP-29 [Fe/H] = 0.11 dex, Hellieretal. 20f 



HD 149026 [Fe/H ] = 0.36 dex, 
[Fe/H] = 0.12 dex. iHolman etal 



Sato et al. 2005, and Kepler-9 
I20T0 ). Figure |6] shows the ob- 



served trend between !??-[Fe/H], and highlights the large diver- 
sity of planetary internal structure and composition. WASP-39's 
metallicity ([Fe/H] = -0.12 dex) strengthens this correlation 
possibly supporting the core-accretion scenario f or planet for- 
mation (iGuillot etal J 120061: iHartman e t al."2009': 'Bouchv et dj 
201(J. The Ohmic dissipation mechanism (Batvgin et al. 201 
coupled to the diversity of heavy element abundances (planetary 
structural evolution), could provide a universal mechanism able 
to slow down planetary cooling and thus the shrinking of the ra- 
dius, possibly explaining the observed correlation between plan- 
etary radii and equilibrium temperature. Further progress may 
be achieved with a better understanding of the structure and at- 
mospheric composition of transiting planets (clouds, chemistry, 
etc). 
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